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I.  INTRODUCTION 


In  a  recent  report.  It  was  shown  that  the  bond-energy  bond-order  (BEBO) 
method  »  *  provided  good  agreement  between  measured  and  computed  rate  coeffi¬ 
cients  for  the  reaction  of  F  with  H2  or  D2.  Additionally,  population  inver¬ 
sions  could  be  predicted  with  this  method  because  the  computed  rate  coeffi¬ 
cients  for  the  production  of  vibrationally  excited  HF  and  DF  were  larger  than 
for  ground-state  HF  and  DF.*  Consequently,  we  investigated  the  applicability 
of  the  BEBO  method  for  estimating  the  effect  of  reactant  vibrational  excita¬ 
tion  on  enhancing  the  kinetics  of  such  reactions  as 


A  +  HB(v  «*  0,1,2,  ...)  =  B  +  AH 


Reactions  of  this  type  have  been  of  considerable  importance  in  connection  with 
de-excitation  in  high-power  chemical  lasers,5’^  isotope  separation  by  selec- 
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tive  laser  excitation,  and  research  in  excited-state  kinetics.  ’» 

2 

The  BEBO  method  was  developed  to  provide  convenient  and  rapid  computa¬ 
tions  for  estimating  rate  coefficients  and  activation  energies  of  hydrogen- 
transfer  reactions.  The  method  utilizes  known  molecular  parameters  such  as 
bond-dissociation  energies,  bond  lengths,  and  vibration  frequencies;  it  is 
based  on  the  principles  of  small-vibration  dynamics  in  molecules  and  on 
transition-state  theory;  and  it  has  provided  satisfactory  agreement  with  many 
measured  rate  coefficients. No  adjustable  parameters  are  required  with 
this  method,  and  it  is  much  briefer  computationally  than  ab  initio  quantum 
mechanical  methods**  or  the  London-Eyring-Polanyi-Sato  (LEPS)  transition-state 

procedure,  which  requires  calculations  for  a  relatively  large  potential  energy 
12 

surface- 


II.  COMPUTATION  METHOD  FOR  REACTIONS  INVOLVING  EXCITED  STATES 


The  BEBO  method  is  based  on  the  concept  that  energy  released  as  the  product 
AH  Is  formed  can  supply  part  or  all  of  the  energy  needed  during  the  breaking  of 
the  reactant  HB  bond.  If  the  energy  available  as  the  bond  order  of  product  AH 
increases  from  zero  to  one  is  less  than  the  energy  needed  as  the  bond  order  of 
reactant  HB  decreases  from  one  to  zero,  the  difference  is  computed  as  the  po¬ 
tential  energy  of  activation.  This  difference  is  calculated  by  the  computer 

2 

program  on  the  basis  of  the  following  BEBO  method  trial  assumptions: 

1.  The  Pauling  rule  for  the  relationship  between  bond  length  and  bond 
order  applies  in  the  transition  state. 

2.  The  Lennard-Jones  6-12  energy  and  distance  parameters  for  the  noble 
gases  represent  the  zero-order  bond  energies  and  lnternuclear  dis¬ 
tances  for  elements  in  the  corresponding  row  of  the  periodic  table. 

3.  The  triplet  repulsion  that  results  from  uncoupled  spins  is  repre¬ 
sented  by  a  Sato  anti-Morse  function  similar  to  that  for  the  triplet 
3  + 

E  state  of  H*. 
u  L 

4.  The  logarithm  of  bond  energy  is  proportional  to  the  logarithm  of  bond 
order. 

The  Sato  triplet  repulsion  is  a  function  of  the  bond  dissociation  energy 
of  the  compound  AB.  For  a  typical  excited-state  reaction  such  as  H  +  HF(v  *  1) 

■  F  +  H2,  note  that  the  bond-dissociation  energy  of  the  vibrationally  excited 
reagent  is  considerably  less  than  that  for  the  ground-state  reagent.  Conse¬ 
quently,  a  bond-dissociation  energy  for  the  triplet  repulsion  that  is  propor¬ 
tionately  lower  than  that  for  the  ground  state  is  used  by  the  BEBO  computations 

for  excited  states.  For  example,  in  'he  computation  for  this  reaction  of  an  H 

atom  with  HF  (v  *  1)  to  produce  ground-state  H2,  the  dissociation  energy  used 
for  triplet  repulsion  is  that  for  HF(v  -  1),  namely,  123.8  kcal/mole,  as  given 

in  Table  l.  In  this  example,  the  value  of  103.3  kcal/mole  would  be  used  as  the 

dissociation  energy  of  H2.  Bond  dissociation  energies  of  excited  states  were 
obtained  from  Reference  6.  Molecular  parameters  used  for  these  BEBO  calcula¬ 
tions  are  listed  in  Table  1. 
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Table  1.  Molecular  Parameters 


Reagent 

Bond 

Dissociation 

Energy,3 

kcal/mol 

Bond  Dissociation 
Energy  for  First 
Vibrational  Level, 8 
kcal/mol 

Interatomic 
Distance,0 
10'8  cm 

Vibrational 

Wave 

Number,0 

cm- 

OH 

101.3 

91.1 

0.971 

3735 

HBr 

86.6 

79.3 

1.41 

2650 

HC1 

102.3 

94.1 

1.28 

2990 

HF 

135.1 

123.8 

0.917 

4138 

DF 

135.7 

127.4 

0.917 

2998 

h2 

103.3 

91.4 

0.742 

4405 

d2 

105.0 

96.5 

0.742 

3119 

aRef.  13. 
^Ref.  6. 
cRef.  14. 


The  trial  assumptions  have  provided  satisfactory  agreement  between  com¬ 
puted  and  experimental  rate  coefficients  and  activation  energies  for  several 
1-4 

reactions.  Nevertheless,  such  assumptions  are  primarily  empirical,  and 
calculations  made  on  the  basis  of  these  assumptions  must  be  considered  to  be 
only  estimates  of  the  rate  coefficients  and  activation  energies  for  the 
reactions  under  consideration. 
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III.  RESULTS  AND  DISCUSSION 


A.  Br  + 

The  rate-enhancing  effect  of  vibrational  excitation  of  H2  to  the  first 
vibrational  level  was  the  subject  of  trajectory  calculations  by  Sims,  Dosser, 

q 

and  Wilson  for  the  reaction 

Br  +  H2(v  -  0,  1)  -  H  +  HBr  (1) 

at  600  and  800  K.  The  results  of  these  calculations  were  3.7  x  10^  at  600  K 
and  3.1  x  10®  at  800  K  for  kj/kg,  the  ratio  of  the  rate  coefficient  for  v  »  1 
to  the  rate  coefficient  for  v  «  0.  The  corresponding  values  of  kj/kg  obtained 
by  our  BEBO  computation  procedure  are  2.9  x  10^  at  600  K  and  2.1  x  10®  at  800 
K.  According  to  the  BEBO  computations,  furthermore,  it  is  predicted  that,  at 
lower  temperatures,  the  rate-enhancing  effect  of  H2  excitation  would  be 
greater.  Thus,  the  computed  BEBO  values  of  kj/Rg  are  9  x  10®,  5  x  10^,  and 
2  x  i05  at  300,  400,  and  500  K,  respectively.  This  trend  results  because  the 
computed  BEBO  activation  energy  decreases  to  8.3  kcal/mole  for  v  =•  1,  compared 
to  20.7  kcal/mole  for  v  «  0. 

B.  Br  +  HC1 

The  kinetics  effect  of  vibrational  excitation  of  HC1  by  a  pulsed  chemical 
laser  was  measured  by  Arnoldi  and  Wolf rum1 ®  for  the  reaction  of  HC1  with  Br 

Br  +  HCl(v  -  0,  1,  2)  -  Cl  +  HBr  (2) 

They  observed  a  vibrational  enhancement*®’*®  of  2  x  10®  for  k^/kg.  The  com¬ 
puted  BEBO  kj/kg  is  1.53  x  10®.  The  measured*’®  enhancement  for  the  second 
vibrational  level  of  HC1  was  1  x  10**  for  k2/kg:  the  BEBO  result  is  5  x  10**. 
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The  computed  BEBO  activation  energy  for  HCl(v  =  2)  in  this  reaction  was 
0.5  kcal/mole,  in  contrast  to  BEBO  activation  energies  of  8.6  kcal/mole  and 
16.9  kcal/mole  for  HCl(v  =1)  and  HCl(v  =  0),  respectively.  It  has  previously 
been  indicated  that  the  BEBO  method  cannot  be  expected  to  provide  reliable 
computations  of  rate  coefficients  for  certain  reactions  in  which  the  bending- 
force  constant  of  the  transition  state  is  so  low  that  the  bending  amplitude  of 
the  transition  state  becomes  unrealistically  high  at  temperatures  above  273  K. 
The  BEBO  calculations  for  HCl(v  =  3,  4,  5)  in  reaction  (2)  revealed  that  this 
bending-force  constant  did  become  too  low  to  permit  reliable  calculations  of 
the  corresponding  rate  coefficients.  This  effect  was  encountered  in  calcula¬ 
tions  for  very  high  vibrational  excitation  in  many  of  the  reactions,  which 
will  be  discussed  in  this  report.  When  the  reagent  vibrational  excitation  is 
sufficiently  high,  the  height  of  the  potential  energy  barrier  approaches  zero; 
consequently,  the  transition  state  calculations  become  unreliable. 

C.  F  +  HC1 

As  part  of  an  extensive  series  of  pioneering  studies  in  Polanyi’s  labora¬ 
tory,  measurements^’^  were  made  to  determine  the  effect  of  HC1  vibrational 
excitation  on  the  reaction 


F  +  HCl(v  -  0,  1)  =  Cl  +  HF 


The  measured  kj/kg  was  in  the  range  of  3  to  9;  the  computed  BEBO  value  of 
kj/kg  is  3.5,  in  sufficient  agreement.  The  BEBO  activation  energies  are  1.4 
and  1.2  kcal/mole  for  v  =  0  or  1.  Reactions  (1)  and  (2)  are  strongly  endo¬ 
thermic  for  v  ■  0,  whereas  reaction  (3)  is  exothermic.  Nevertheless,  the 
computed  BEBO  results  are  generally  in  good  agreement  with  the  measured  or 
trajectory  values  for  rate  enhancement  owing  to  vibrational  excitation.  It 
has  previously  *  been  explained  that  the  BEBO  method  is  based  on  a  model 


BH  +  A 


B  +  HA 


i 


where,  in  the  transition  state  B  ...  H  ...  A,  the  bond-order  m  of  the  new  bond 
is  increasing  from  zero  to  one,  while  the  bond  n  of  the  breaking  bond  B...H  is 
simultaneously  decreasing  from  one  to  zero.  With  the  BEBO  model,  the  energy 
available  from  the  formation  of  the  H...A  bond  can  be  used  to  supply  part  or 
all  of  the  energy  required  for  dissociating  the  B...H  bond.  A  potential- 
energy  barrier  develops  for  most  H-transfer  reactions  because  the  energy 
available  during  the  formation  of  H...A  does  not  supply  all  the  energy  needed 
during  the  dissociation  process  for  the  reagent  B...H.  However,  this  BEBO 
model  suggests  that,  if  the  reagent  is  vibrationally  excited  and  consequently 
requires  less  energy  for  its  dissociation,  the  resultant  activation  energy 
barrier  is  expected  to  be  lower  than  for  the  case  when  v  *  0  for  reagent  BH. 
The  BEBO  computations  were  generally  in  agreement  with  this  expectation.  The 
explanation  of  many  cases  of  enhanced  rates  resulting  from  reagent  vibrational 
excitation  would,  therefore,  simply  be  based  on  the  lower  bond  dissociation 
energy  of  the  excited  reagent  plus  the  coordinated  availability  of  energy  from 
the  formation  of  the  new  bond  H...A.  In  the  BEBO  model,  this  coordination 
arises  from  maintaining  the  sum  of  the  bond  orders  n  and  m  equal  to  1.00, 
inasmuch  as  one  bonding  orbital  on  H  is  transferring  from  B  to  A  during  the 
reaction. 

0.  H  +  H0 

Heidner  and  Kasper*®  measured  the  predominantly  reactive  Increase  in  rate 
coefficient  resulting  from  excitation  of  H2  for 

H  +  H2(v  -  0,  1)  -  H  +  H2  (5) 

They  obtained  a  value  of  about  1500  for  k^/kp  at  300  K.  The  computed  BEBO 
value  for  k^/kp  is  6000,  and  the  activation  energy  for  v  ■  1  is  3.8  kcal/mole 
compared  to  9.0  kcal/mole  for  v  ■  0. 


13 


For  the  reaction 


0  +  H2(v  =0,  1)  =  H  +  OH  (6) 

Light^^  measured  k^/kp  to  be  2600.  The  computed  BEBO  k^/kg  is  5400,  and  the 
computed  activation  energies  are  5.7  and  10.8  kcal/mole  for  v  =*  1  and  0, 
respectively. 

F.  0  +  HC1 

20 

MacDonald  and  Moore  measured  the  rate  constant  for  the  reaction 

0  +  HCl(v  =  1)  =  Cl  +  OH  (7) 

and  found  it  to  be  5  x  10^  citrV mole-sec.  The  computed  BEBO  rate  constant  is 
3  *  1011  cm^ /mole-sec.  The  computed  BEBO  value  of  k^/kQ  is  45,  with  the 
energies  of  activation  equal  to  2.6  and  1.8  kcal/mole  for  v  -  0  and  1, 
respectively. 

G.  D  +  HC1 

21 

Kneba  used  time-resolved  mass  spectroscopy  and  atomic  resonance  absorp¬ 
tion  to  determine  the  rate  coefficient  at  298  K  for  the  reaction  of  HCl(v  *  1) 
with  D  to  produce  DH  and  Cl.  The  measured  rate  coefficient  was  2.6  x  10^ 
cm^/mole-sec,  whereas  the  computed  BEBO  rate  coefficient  at  300  K  is  2  x  10*^ 
cnr /mole -sec.  The  BEBO  activation  energy  is  3.0  kcal/mole. 

H.  I  +  HC1 

Kaplan,  Levine,  and  Manz^  applied  their  information  theory  (surprisal) 
method  to  computing  rate  coefficients  at  300  K  for  the  reaction  of  HCl(v)  with 
and  iodine  atom  to  form  HI  and  Cl,  making  use  of  the  experimental  rate  coeffi¬ 
cients  for  the  reverse  reactions.  The  computed  BEBO  rate  coefficients  are  in 
fair  agreement  with  their  results.  For  values  of  v  ■  0,  1,  2,  3,  and  4,  the 
ratios  of  the  BEBO  rate  coefficients  at  300  K  to  their  surprisal  rate  coeffi¬ 
cients  are,  respectively,  2.9,  3.2,  1.3,  1.1,  and  2.2.  Figure  1  reveals  the 
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Figure  1.  Dependence  of  log  of  BEBO  Rate  Coefficient  on  Reagent  Excitation 


I 


dependence  of  the  BEBO  rate  coefficients  on  the  vibrational  excitation  of 
HC1.  There  is  an  almost  linear  dependence  of  log  ky  on  the  vibrational  exci¬ 
tation  of  HCl(v)  up  to  a  v  of  4.  This  approximate  linearity  occurs  because 
the  computed  BEBO  decrease  in  the  height  of  the  potential  energy  of  activation 
barrier  for  v  »  1,  2,  3,  or  4  corresponds  very  closely  (average  deviation  only 
0.11  kcal/mole)  with  the  increase  in  HCl  vibrational  excitation  energy  (Cohen 
and  Bott^)  for  v  *  1,  2,  3,  and  4.  The  reaction  becomes  exoergic,  however, 
for  v  -  5,  and  the  BEBO  decrease  in  barrier  height  is  typically  much  less  than 
the  increase  in  reagent  excitation  for  exoergic  reactions  for  HCl(v  *  0,  1,  2, 
3,  4). 

The  importance  of  minimizing  reactive  de-excitation  of  vibrationally  ex¬ 
cited  HF  or  DF  is  well  recognized^  by  workers  concerned  with  increasing  the 
power  from  HF  or  DF  laser.  The  BEBO  method  can  provide  a  convenient  and  eco¬ 
nomical  way  to  estimate  rate  coefficients  for  such  reactions.  When  measured 
de-excitation  rate  coefficients  do  not  exist  for  the  de-excitations  of  inte¬ 
rest,  the  BEBO  results  can  be  utilized  in  the  several  available  computer  codes 
for  modeling  the  behavior  of  chemical  lasers.  The  results  of  a  typical  BEBO 
de-excitation  computation  are  given  in  Table  2  for  the  chemical  laser  reaction 
in  which  excited  HF  is  lost  by  reaction  with  H 


H  +  HF(v  =  0  -  10)  =*  F  +  H2(v'  «  0,1) 


The  BEBO  method  uses  so  little  computer  time  that  it  is  economical  to  make 
calculations  for  values  of  v  as  high  as  10  and  also  to  include  reactions  in 
which  the  product  is  excited  to  v'  =  1.  Furthermore,  it  was  a  standard  pro¬ 
cedure  in  this  series  of  computations  to  obtain  the  results  for  temperatures 
of  100,  400,  500,  600,  800,  1000,  1200,  and  1400  K,  as  well  as  300  K.  Values 
were  also  computed  for  the  activation  energy  (Ea)  and  the  potential  energy  of 
activation  (V*).  In  Table  2,  AB  is  the  bond  dissociation  energy  of  the  re¬ 
agent  minus  the  bond  dissociation  energy  of  the  product. 
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Table  2.  BEBO  Estimates  of  Bate  Coefficients  for 
H  +  HF(v  -  0  -  10)  «  F  +  H2(v'  -  0,  1) 


Rate  Coefficients  at  300  K, 
cm3/mol-sec 


V 

V 

-  0 

V 

-  1 

kcal^mole 

V,* 

kca 1/mole 

AB, 

kcal/mole 

0 

0.41 

X 

10" 11 

0.38 

X 

10-19 

34.6 

33.2 

+31.8 

1 

0.91 

X 

10~3 

0.11 

X 

10-10 

23.0 

21.9 

+20.5 

2 

0.11 

X 

106 

0.13 

X 

IM 

1 

o 

11.7 

11.8 

+  9.6 

3 

0.40 

X 

1012 

0.48 

X 

105 

2.7 

4.4 

-  0.7 

4 

0.11 

X 

1014 

0.21 

X 

1013 

1.3 

1.1 

-10.6 

5 

0.78 

X 

1014 

0.30 

X 

1014 

1.2 

0.2 

-20.0 

6 

0.73 

X 

1015a 

0.33 

X 

1015a 

7 

0.34 

X 

1016a 

0.34 

X 

1016a 

8 

0.61 

X 

1016a 

0.61 

X 

l016a 

9 

0.12 

X 

1017a 

0.11 

X 

10l7a 

10 

0.26 

X 

1017a 

0.25 

X 

l017a 

aAt  these  values  of  v,  the  computed  values  of  the  bending-force  constant  in 
the  transition  state  were  so  low  that  small  vibration  dynamics  (and  BEBO) 
do  not  apply.2*® 

The  results  given  in  Table  2  are  typical  of  reactions  that  are  strongly 
endothermic  when  the  reagent  is  in  the  ground  state.  For  v  »  0,  the  corre¬ 
sponding  rate  coefficient  (kp)  is  so  small  that  rates  would  be  extremely  slow. 
For  v  ■  1,  the  rate  coefficient  (kj)  increases  by  a  factor  of  more  than 
2  *  10®,  but  the  rate  of  de-excitation  of  HF(v  *  1)  by  reaction  with  H  would 
be  small  on  the  time  scale  of  a  typical  chemical  laser.  However,  for  HF  ex¬ 
cited  to  v  «  4  or  5,  the  rate  coefficients  for  de-excitation  by  reaction  with 
H  are  quite  large  and  are  significant  for  affecting  the  behavior  of  HF  lasers. 
Depending  on  the  characteristics  of  the  chemical  laser,  the  de-excitation  of 
HF(v  -  3)  by  H  might  be  appreciable.  In  Table  2,  it  is  evident  that  AB,  Efl, 
and  V*  decrease  when  v  Increases.  As  the  reaction  becomes  more  exothermic 
(for  example,  at  v  *  4),  the  changes  in  rate  coefficient  and  activation  energy 
become  smaller.  In  Fig.  1,  the  dependence  of  log  by  on  HF  excitation  is 


17 


fairly  linear  for  v  ■  1  or  2.  However,  when  AB  becomes  negative  at  higher 

values  of  v,  the  Increases  In  log  kv  become  comparatively  smaller  because  the 

decreases  In  BEBO  Ea  and  V*  are  less  than  the  corresponding  decreases  In  AB. 

For  v  *  4  or  5  (Table  2),  the  rate  coefficients  for  the  formation  of  vibra- 

tionally  excited  product  H2(v  =  1)  may  be  significant  also.  In  Table  2,  the 

superscript  a  is  used  to  denote  those  cases  (at  relatively  high  values  of  v 

and  exothermicity)  for  which  the  computed  bending  force  constants  in  the 

transition  state  are  so  low  that  the  amplitudes  of  vibration  in  the  transition 

2  3 

state  are  unrealistically  *  large  and  the  potential  energy  of  activation  is 

virtually  nil.  In  such  cases  BEBO  rate  coefficients  should  not  be  used  be- 

2 

cause  BEBO  relies  on  transition  state  theory;  Instead,  it  is  believed  prefer¬ 
able  to  use  collision  theory  to  calculate  rate  coefficients. 

Table  3  is  a  summary  of  the  BEBO  results  for  computations  similar  to 
those  made  for  Table  2.  Any  oxygen  atoms  present  in  the  chemical  laser  gases 
would  have  HF  de-excitation  rate  coefficients  somewhat  larger  than  those  for 
hydrogen  atoms  at  v  higher  than  2.  An  atom,  such  as  Br,  which  reacts  more 
endothermically  with  HF(v),  would  generate  lower  de-excitation  rate  coeffi¬ 
cients.  When  v  is  5  or  higher,  however,  it  appears  that  almost  any  atom  or 
radical  present  in  the  laser  gas  could  react  rapidly  with  HF(v)  to  de-excite 
it.  Bernstein  has  observed1®  that  reactions  that  are 


Table  3.  Rate  Coefficients3  for  De-excitation  of  HF(v)  and  HCl(v) 


V 

0  +  HF(v) 

-*•  F  +  OH 

Br  +  HF(v) 

-*■  F  +  HBr 

H  +  HCl(v) 

>  Cl  +  h2 

0 

1.1  x 

10" 12 

5.2  x 

10"?A 

8.5 

x  109 

i 

2.4  x 

10I4 

1.4  x 

10" 14 

2.3 

x  1011 

2 

2.7  x 

104 

1.9  x 

10~6 

2.2 

x  10  2 

3 

6.4  x 

10 \\ 

9.2  x 

l0Q 

1.3 

x  1013 

4 

2.2  x 

1013 

2.2  x 

l0?/ 

4.4 

x  1013 

5 

1.3  x 

1014 

5.0  x 

10  4b 

1.5 

x  1014 

6 

1.7  x 

1015b 

8.0  x 

10l4b 

7.1 

x  1014b 

aRate  coefficients  at  300  K,  cm3/mole-sec. 

^Bending  force  constant  of  transition  state  is  too  low. 
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strongly  endothermic  for  ground-state  HF  have  a  vibrational  threshold  effect. 
The  BEBO  results  for  HF  are  consistent  with  this  observation  because  the  rate 
coefficients  for  the  lower  values  of  v  In  Tables  2  and  3  are  so  small  that  the 
experimentally  measured  rates  could  be  extremely  low  compared  to  the  rates  for 
v  =  4,  5,  or  6.  For  the  more  thermoneutral  reaction  of  H  with  HC1,  however. 

It  Is  evident  In  Table  3  that  the  rate  coefficients  do  not  Increase  so  rapidly 
as  a  function  of  v. 

In  Table  4,  results  are  summarized  for  BEBO  computations  involving  DF(v). 
The  removal  rates  for  HF(v)  are  important,  of  course,  for  DF  lasers.**  A  com¬ 
parison  of  the  rate  coefficients  In  Table  2  with  those  for  H  +  DF(v)  (Table  4) 
at  equal  values  of  v  indicates  that  the  removal  rate  coefficients  for  H  +  DF 
are  generally  lower  than  those  for  H  +  HF(v).  This  difference  Is  probably 
largely  the  result  of  the  lower  exothermlcity  for  the  H  +  DF(v)  reactions  be¬ 
cause  of  the  closer  spacing  of  the  DF  vibrational  levels.  The  larger  computed 
BEBO  rate  coefficients  for  the  D  +  DF(v)  reactions  can  be  attributed  to  the 
lower  vibration  frequency  of  D2  compared  to  that  for  DH.  The  calculated  de- 
excltatlon  rate  coefficients  for  D  +  DF(v)  at  v  -  3,  4,  5,  or  6  are 
considerably  smaller  than  those  for  H  +  HF(v).  This  factor  could  contribute 
to  higher  laser  power  for  the  DF  laser  than  for  the  HF  chemical  laser  at 
higher  vibrational  transitions. 

In  the  last  two  columns  of  Table  4  are  the  rate  coefficients  for  DF(v) 
de-excitation  by  means  of  exchange  reactions  involving  D  atom  transfer.  Note 
that  the  rate  coefficients  for  de-excitatlon  of  DF(v)  by  F  are  considerably 
higher  than  those  for  DF(v)  de-excltation  by  D.  This  fact  might  indicate  that 
obtaining  higher  power  from  DF  chemical  lasers  would  be  facilitated  by  mini¬ 
mizing  the  quantity  of  F-atoms  present  in  the  laser  gases  after  the  F  +  D2  ■ 
DF(v)  +  F  results  in  rapid  de-excltatlon  by  the  nearby  F  atoms  exchanging  with 
the  DF(v).  The  high  endothermicity  of  the  D  +  DF  reaction  is  the  probable 
reason  that  these  rate  coefficients  are  smaller  than  those  for  the  F  +  DF 
reaction.  The  BEBO  results  apply  only  to  reactive  and  not  to  nonreactive  de- 
excitatlons.  The  values  In  the  last  column  of  Table  4  differ  from  those  of 
the  next-to-last  column  because  the  DF  of  the  last  column  has  been  de-excited 
to  the  v  -  1  level  rather  than  to  the  v  ■  0  level.  In  the  BEBO  method,  it  is 
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Table  4.  De-excitatlon  Rate  Coefficients3  for  DF(v) 


D  +  DF(v)  H  +  DF(  v)  F  +  DF(v)  F  +  DF(v) 

v  r  F  +  D2  ~  F  +  HD  2:  FD(v  -  0)  +  F  ~  FD(v  -  1)  +  F 


0 

0.38 

X 

►— 

o 

o 

0.39 

X 

10-11 

1 

0.46 

X 

10"4 

0.58 

X 

10"5 

2 

0.34 

X 

102 

0.42 

X 

10 

3 

0.12 

X 

00 

o 

0.16 

X 

107 

4 

0.17 

X 

1012 

0.36 

X 

1011 

5 

0.56 

X 

1013 

0.27 

X 

1013 

6 

0.29 

X 

1014 

0.18 

X 

1014 

7 

0.13 

X 

10« 

0.82 

X 

1014 

8 

0.10 

X 

1016b 

0.50 

X 

1015b 

9 

0.42 

X 

l016b 

0.27 

X 

1016b 

0.13 

X 

1012 

- 

— 

0.44 

X 

1013 

0.44 

X 

1013 

0.19 

X 

1014 

0.10 

X 

1014 

0.17 

X 

1015 

0.19 

X 

1014 

0.16 

X 

10l6b 

0.77 

X 

1014 

0.22 

X 

1016b 

0.22 

X 

l015b 

0.26 

X 

1016b 

0.30 

X 

1015b 

3Rate  coefficients  at  300  K,  cra3/mole-sec . 
bBending-force  constant  of  transition  state  Is  too  low. 


to  be  expected  that  the  rate  coefficients  for  DF(v  -  1)  production  would  be 
smaller  than  those  for  ground-state  DF  production  because  the  bond  energy  is 
smaller  for  DF(v  -  1)  than  for  ground  state  DF  If  v  is  at  least  2.  However, 
consideration  of  the  transition  state  F...D...F  for  this  reaction  Indicates 
that  extra  energy  may  be  required  to  dissociate  FDF  because  of  the  formation 
of  a  hydrogen  bond  such  as  that  occurring  in  hydrogen  fluoride  vapor. 23  This 
extra  bond  energy  can  range  from  2  to  10  kcal/mole.  BEBO  calculations  were 
made  In  which  this  possible  additional  stability  for  FDF  was  considered;  the 
calculated  rate  coefficients  were  found  to  be  lowered  by  a  factor  of  one 
hundred  for  10  kcal/mole  stabilization.  Such  extra  bond  energy  for  FDF  could 
possibly  be  the  reason  that  the  BEBO  rate  coefficients  are  higher  (Table  4) 
for  the  loss  of  several  vibrational  quanta  than  for  th  loss  of  one  quantum, 
whereas  calculations  by  Wilkins24  indicate  higher  rate  coefficients  when  only 
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one  vibrational  quantum  is  lost.  BEBO  computations  were  also  made  for  de- 
excitation  of  HF(v)  by  exchange  reactions.  For  the  reaction  F  +  HF(v  ■  1)  *  F 
+  HF,  the  calculated  rate  coefficient  was  0.13  *  10^  at  300  K.  This  value  is 
much  higher  than  for  HF(v  *  1)  de-excitation  by  H  (Table  2).  For  v  >  1,  the 
bending-force  constants  were  too  low;  this  Indicates  much  faster  de-excitation 
of  HF(v)  by  F  atom  exchange  than  by  reaction  with  H  atoms. 
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IV.  BEBO  CALCULATIONS  FOR  DEUTERIUM  ISOTOPE  SEPARATION  REACTIONS 
INVOLVING  VIBRATIONALLY  EXCITED  REAGENTS 


Because  of  the  desirability  of  developing  more  economical  methods  of 
producing  D2O  for  use  in  heavy-water  moderated  nuclear  power  reactors,  there 
has  been  speculation  on  the  bimolecular  reaction  rate  enhancement  obtainable 
by  selective  excitation  of  deuterium  compounds.2^  In  Table  4,  note  that  the 
BEBO  method  estimate  for  the  rate  coefficient  of  the  reaction  of  DF(v  -  5)  with 
H  to  produce  DH  is  0.27  *  10^  cm^/mole-sec,  whereas  the  rate  coefficient  of 
the  same  reaction  for  ground  state  DF  or  HF  is  less  than  0.5  x  10_1J  cm^ /mole- 
sec.  Because  of  the  endothermicity  of  this  reaction,  however,  it  would  be  ne¬ 
cessary  to  expend  enough  energy  to  excite  DF  to  its  fifth  vibrational  level  in 
order  to  achieve  this  desirable  enhancement  of  the  rate  coefficient.  BEBO 
method  results  for  less  endothermic  reactions  of  deuterium  compounds  are  sum¬ 
marized  in  Table  5.  These  reactions  do  not  require  as  much  excitation  of  the 
deuterium  compounds  as  do  more  endothermic  reactions,  in  order  to  achieve  rates 
that  could  yield  adequate  productivity.  As  an  example  (Table  5),  DF(v  *  2) 
reacting  with  OH  has  a  rate  coefficient  of  0.22  *  10^2  cm^/mole-sec,  compared  to 
4.2  cm  /mole-sec  for  the  H  +  DF(v  -  2)  reaction.  The  BEBO  rate  coefficients  for 
OH  +  DF  with  v  equal  to  3  or  4  were  0.30  *  10*^  or  0.33  x  10^  cm^/mole-sec, 
respectively.  These  k^/s  are  also  much  higher  than  the  corresponding  k^/s 
for  H  +  DF  in  Table  4.  The  decrease  in  &B  (the  difference  in  bond  dissocia¬ 
tion  energies  for  product  and  reactant)  for  the  OH  +  DF  reaction  is  16.4 
kcal/mole  when  v  changes  from  0  to  2  (Table  5).  The  corresponding  decrease  in 
activation  energy  is  16.3  kcal/mole  as  computed  by  the  BEBO  method.  This 
computation  reveals  that  the  vibrational  excitation  of  DF  can  be  used  effi¬ 
ciently  to  lower  the  activation  energy  barrier  for  DF  reactions.  However,  at 
higher  values  of  v  the  decrease  in  reaction  endothermicity  may  reduce  the 
efficiency  of  vibrational  excitation  for  lowering  the  activation  energy.  For 
example,  in  the  OH  +  DF  reaction,  Ab  decreases  by  15.3  kcal/mole  when  v  rises 
from  2  to  4,  but  the  computed  E_  drops  by  only  1.7  kcal/mole. 
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Table  5.  BEBO  Rate  Coefficients3  for  Reactions  of  Vibrationally  Excited  Deuterium  Compounds 


Because  the  bond  dissociation  energy  of  DC1  is  less  than  that  for  DF, 
higher  rate  coefficients  (and  greater  rates  of  formation  of  deuterium  compound 
products)  could  be  obtained  at  equal  values  of  v  although  less  energy  would  be 
needed  to  excite  DC1  than  DF.  The  BEBO  results  in  Table  5  for  the  0  +  DC1  and 
H  +  DC1  reactions  support  this  possibility.  Thus,  DF  must  be  excited  to 
v  4  to  obtain  (Table  4)  a  BEBO  kr  of  0.36  x  10^  for  the  H  +  DF  reaction, 
but  a  BEBO  kr  of  0.87  x  10^  is  obtained  at  v  ■  1  for  the  H  +  DC1  reaction 
(Table  5).  Vibrational  excitation  of  DC1  to  v  *  1  appears  to  be  quite  effi¬ 
cient  for  lowering  the  activation  energy  barrier  of  the  0  +  DC1  reaction  be¬ 
cause  the  BEBO  Ea  is  decreased  by  6.1  kcal/mole  (Table  5)  as  a  result  of 
DCl(v  *  0)  being  excited  by  6.0  kcal/mole  to  v  *  1.  It  is  possible  that  a  DC1 
discharge  lamp  could  be  a  satisfactory  source  of  radiation  for  selectively 
exciting  DC1  (in  a  natural  mixture  of  DC1  and  HC1)  to  v  *  1  with  rotational 
quantum  numbers  from  0  to  10.  The  BEBO  k  's  indicate  that  this  excitation 

^  y  r 

would  provide  rate  coefficients  for  the  D  +  DC1  reaction  that  are  10^  to  10° 
larger  than  the  k^'s  for  reaction  with  ground-state  DC1  or  HC1. 

Inasmuch  as  the  H  +  DC1  reaction  is  exothermic,  vibrational  excitation 
does  not  cause  as  large  an  Increase  in  the  BEBO  kr's  as  it  does  for  the  0  +  DC1 
reaction.  Vibrational  excitation  of  the  D-C  bond  in  DCFj  or  DCCl^  can  also  be 
effective  for  enhancing  rate  constants,  as  can  be  seen  in  the  BEBO  results  of 
Table  5.  For  Cl  +  DCF3,  the  computed  activation  energy  barrier  was  lowered  by 
4.5  kcal/mole  when  the  D-C  bond  was  excited  by  4.6  kcal/mole  to  v  *  I.  The 
BEBO  k^’s  for  0  +  DCF3  indicate  that  this  reaction  is  too  slow  for  satisfac¬ 
tory  productivity  of  OD  even  when  v  »  2  for  the  D-C  bond.  On  the  other  hand, 
the  BEBO  kr  for  the  Br  +  DCCI3  reaction  is  rather  large  even  when  v  «  1  for 
the  D-C  bond.  Also,  the  computed  Efl  for  this  reaction  decreases  by  3.3 
kcal/mole  when  the  D-C  bond  is  excited  by  4.2  kcal/mole  to  v  *  1. 


25/o?6 


1 


REFERENCES  ,  /  /' 

rr 

r 

1.  S.  W.  Mayer  and  N.  Cohen,  Rate  Coefficients  Calculated  by  the  Bond-Energy 
Bond-Order  Method  for  Excited  State  F  +  H?  °  H  +  HF  Reactions, 
TR-0079(4940-01)-2,  The  Aerospace  Corp.,  El  Segundo,  Calif.  (16  February 
1979). 

2.  H.  S.  Johnston,  Gas  Phase  Reaction  Rate  Theory,  The  Ronald  Press  Co.,  New 
York  (1966). 

3.  S.  W.  Mayer  and  L.  Schieler,  J.  Phys.  Chem.  72,  236  (1968). 

4.  S.  W.  Mayer,  J.  Phys.  Chem.  73,  3941  (1969). 

5.  D.  J.  Spencer,  T.  A.  Jacobs,  H.  Mirels,  and  R.  W.  F.  Gross,  Int.  J.  Chem. 
Klnet .  J j  493  (1969);  also,  Appl.  Phys.  Lett.  16,  235  (1970). 

6.  N.  Cohen  and  J.  F.  Bott,  "Kinetics  of  Hydrogen-Halide  Chemical  Lasers," 

Ch.  2  In  R.  W.  F.  Gross  and  J.  F.  Bott,  eds.,  Handbook  of  Chemical 
Lasers,  John  Wiley  and  Sons,  New  York  (1976). 

7.  J.  H.  Birely,  J.  Photochem.  4_,  269-280  (1975). 

8.  J.  C.  Polanyi,  Acc.  Chem.  Res.  5,  161  (1972). 

9.  L.  B.  Sims,  L.  R.  Dosser,  and  P.  S.  Wilson,  Chem.  Phys.  Lett.  32,  150-152 

(1975). 

10.  R.  B.  Bernstein,  "State-to-State  Cross  Sections  and  Rate  Constants  for 
Reactions  of  Neutral  Molecules,"  Ch.  1  in  P.  R.  Brooks  and  E.  F.  Hayes, 
eds.,  State-to-State  Chemistry,  American  Chemical  Society,  Washington, 

D.C.  (1977). 


11.  S.  F.  Boys  and  I.  Shavitt,  J.  Chem.  Phys.  31,  1359  (1959). 

12.  R.  N.  Porter  and  M.  Karplus,  J.  Chem.  Phys.  40,  1105  (1964). 

13.  B.  Darwent,  Bond  Dissociation  Energies  in  Simple  Molecules,  NSRDS-NBS  31, 

National  Bureau  of  Standards,  Washington,  D.C.  (1970). 

14.  D.  R.  Stull  and  H.  Prophet,  JANAF  Thermochemical  Tables,  2nd  ed.,  NSRDS- 
NBS  37,  National  Bureau  of  Standards,  Washington,  D.C.  (1971). 

15.  D.  Arnold!  and  J.  Wolfrum,  Ber.  Bunsenges.  Phys.  Chem.  80,  892  (1976). 

16.  L.  J.  Kirsch  and  J.  C.  Polanyi,  J.  Chem.  Phys.  57,  4498  (1972).  i 

17.  A.  M.  Ding  et  al..  Disc.  Faraday  Soc.  55,  252  (1973).  j 


27 


2mm 


3 


18.  R.  F.  Heidner  III  and  J.  V.  V.  Kasper,  Chem.  Phys.  Lett.  15,  179  (1972). 

19.  G.  C.  Light,  J.  Chem.  Phys.  68,  2831  (1978). 

20.  R.  G.  MacDonald  and  C.  B.  Moore,  J.  Chem.  Phys.  68,  513  (1978). 

21.  M.  Kneba,  Thesis,  Universitat  Gottingen,  1977. 

22.  H.  Kaplan,  R.  D.  Levine,  and  J.  Manz,  Chem .  Phys .  12,  447  (1976). 

23.  L.  Pauling,  The  Nature  of  the  Chemical  Bond,  Chap.  12,  Cornell  University 
Press,  Ithica,  N.  Y.  (1960). 

24.  R.  L.  Wilkins,  J.  Chem.  Phys.  58,  3038  (1973). 

25.  V.  S.  Letokhov,  "Laser  Separation  of  Isotopes,"  in  B.  S.  Rabinovitch, 

ed..  Annual  Review  of  Physical  Chemistry,  Vol.  23,  Annual  Reviews  Inc., 
Palo  Alto,  Calif.  <1977). 


28 


LABORATORY  OPERATIONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver¬ 
satility  and  flexibility  have  been  developed  to  a  high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  apace  and  missile  systems.  Expertise  in  the  latest  scientific  devel¬ 
opments  is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory;  Launch  and  reentry  aerodynamics,  heat  trans¬ 
fer,  reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo¬ 
sensitive  materials  and  sensors,  high  precision  laser  ranging,  and  the  appli¬ 
cation  of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi¬ 
conducting,  superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory;  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor¬ 
rosion  and  fatigue -induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia¬ 
tion  from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia¬ 
tions  in  space  on  space  systems. 
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